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When time reversal symmetry is present there is a sharp distinction between topological and 
trivial band insulators which ensures that, as parameters are varied, these phases are separated by 
a phase transition at which the bulk gap closes. Surprisingly we find that even in the absence of 
time reversal symmetry, gap less regions originating from the phase boundaries persist. Moreover 
the critical line generically opens up to enclose Chern insulating phases that are thin but of finite 
extent in the phase diagram. We explain the topological origin of this effect in terms of quantized 
charge pumping, showing in particular that it is robust to the effect of disorder and interactions. 
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I. INTRODUCTION 

The distinction between Topological and ordinary 
band insulators relies on time reversal symmetry 
(TRS) 1 ' 2 . In the presence of this symmetry the two 
phases must be separated in parameter space by a quan- 
tum phase transition at which the energy gap vanishes. 
Similarly a real space boundary between the two insula- 
tors must sustain protected metallic edge states, as has 
been clearly observed in experiment 3 ' 4 . When time re- 
versal symmetry (TRS) is broken the edge states gen- 
erally become gapped or localized by disorder 5 , albeit 
some residual consequences of the bulk topology, such as 
a surface Chern number, may remain. Because there is at 
least a superficial relation between the edge states in real 
space and the gapless (critical) phase boundaries in pa- 
rameter space, it is natural to ask what is the fate of the 
latter when time reversal symmetry is broken. Are these 
critical regions in the phase diagram immediately gapped 
out or localized, or do they contain some robust features 
that are not destroyed by the symmetry breaking? 

In this paper we asses the robustness of the critical 
lines that enclose a 2D topological insulator in parameter 
space, separating it from an ordinary band insulator. We 
first derive general criterions for stability of those critical 
regions to weak TRS breaking perturbations. We fur- 
ther show that in certain cases the critical regions have 
a stronger topological protection, which does not rely 
on time reversal symmetry and is hence much stronger 
than the protection of the gapped insulators. In this case 
the gapless surface that enclosed the topological phase 
in parameter space may split up into segments to allow 
adiabatic connection between the former topological and 
ordinary insulators. However, the critical segments per- 
sist and cannot terminate for arbitrarily strong symmetry 
breaking perturbations. 

For a two dimensional non interacting and transla- 
tional invariant model, the topological criterion can be 
simply formulated in terms of a winding number on the 
Brillouin zone. The critical point between a topological 
and ordinary band insulator is where the bulk Dirac dis- 
persion becomes massless. These gapless points follow a 
certain trajectory in /c-space as the parameters are varied 



along a critical line that encloses the topological insulator 
in the phase diagram. We assert that the critical regions 
persist for arbitrary time reversal symmetry breaking, if 
the trajectory of Dirac points winds around the Brillouin 
zone. 

To establish the robustness of the critical region also in 
presence of disorder and interactions we make a connec- 
tion between the winding number on the Brillouin zone 
and a Chern number associated with charge pumping. 
Specifically the winding of the Dirac points on the Bril- 
louin zone is directly related to a Chern number defined 
on a path in parameter space that encircles the topologi- 
cal phase from the outside, i.e. on a path that lies entirely 
in the ordinary insulating phase. A non zero Chern num- 
ber implies pumping of an integer number of particles per 
row of the lattice upon following such a closed adiabatic 
path. Hence, the critical phase boundary might exhibit 
a topological protection that is much stronger than that 
of the insulating phases 6 . In contrast to the latter, this 
topological protection does not rely on time reversal sym- 
metry but only on particle number conservation. 



II. ROBUST CRITICAL REGIONS IN AN 
EXTENDED KANE AND MELE MODEL 



To elucidate the main results of the paper we first 
demonstrate the aforementioned features on a concrete 
2D tight binding model which can be tuned between a 
topological insulator (TI) and a band insulator (BI). This 
model is an extension of the Kane and Mele model 1 , de- 
scribing electrons hopping on a Honeycomb lattice in the 
presence of spin-orbit coupling. It includes two param- 
eters, staggered hopping and staggered potential, which 
drive the system into a BI phase. As a function of these, 
the phase diagram consists of an elliptical TI region en- 
circled by a BI phase with a critical line in between. We 
shall study the fate of this line as TRS is broken by a 
staggered and a uniform in-plane Zeeman field. 
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The corresponding Hamiltonian is given by 

H km = H g + V so + V stag + V so + Vsu (1) 

V sta = A^2ticl a Ci a , 

is 

+ i\R^24 s (a ss > x dij) z c js >. 

(ij) 

V 5t =St 4 a (l + sX R /t)c ja 

i-j=±d 3 

~ 6t C ia( 1 - sX R/ t ) C JS, 
i-j=±d 2 

Here q s annihilates an electron with S z = s on the i'th 
site of an hexagonal lattice, & = ±1 is an alternating 
sign on the two sublattices, {{ij)) denotes next-nearest- 
neighbor sites and Vij = 1(— 1) if the lattice path be- 
tween the two sites i and j, winds clock- wise (counter 
clock-wise). The vectors di and the staggered hopping 
pattern are shown in Fig. (1). Summation of repeated 
indices is implicit. Notably the model includes an ex- 
tra Rashba term proportional to St, that couples to the 
out-of-plane spin component. The reason for adding the 
additional spin orbit term is to break spurious symme- 
tries 7 of the model and thus explore generic features. 

In the absence of TRS breaking perturbations, the 
phase diagram of the above model in the plane of (A, St) 
consists of an inner TI phase separated by a critical line 
from an outer BI phase, as shown in Fig. (2a), dashed 
red line. To test the robustness of this line we consider a 
staggered (M stg ) and uniform (M x ) Zeeman field terms 
in the x— direction 

V M = MstgYj&laWxUfT'Ci,*' (2) 

i 

+ M x ^cl a [a x ] a ^Ci^ 

The staggered term may arise if the system breaks 
TRS spontaneously by developing an in-plane anti- 
ferromagnetic order. We have verified that strong 
Hubbard-type repulsive interactions indeed lead to such 
a symmetry breaking with an in-plane magnetization (see 
also 8 ). 

Figure (2) shows the numerically obtained phase di- 
agram in (A, St) plane for different TRS breaking per- 
turbations and fixed values of the remaining parameters 
(t = l,A so = 1 and = 0.5). Fig. (2 a) and (2 
b), show the phase diagram for non-zero staggered field 
(M stg — 2) and non-zero uniform field (M x = 0.5), re- 
spectively. The phase boundary between the TI and BI 
is marked in dashed red line. Interestingly, it does not 
gap out after one breaks TRS. Instead it expands into 



a double line (blue contour) with an IQHE phase with 
(j xy = ±x appearing within. Metallic regions (grey re- 
gions) may also appear due to an indirect closing of the 
gap. 

The gapless band touching lines in Fig. (2) do not 
disappear even as we further increase the TRS breaking 
terms. For M stg = 5t, lOt the eight-shaped contours in 
(2a) simply intersect the A axis further away (at A « 
±M stg ) and gradually become squeezed along this axis. 
For large M x , the closed contour in (2b) eventually splits 
along the St axis and also intersect the A axis near ±M X . 
In addition some new band touching lines appear near 
the original line. This apparent robustness of the band 
touching lines will be explained in the next two sections. 




FIG. 1: Illustration of the lattice model and the different 
couplings given in Eq. (1) 



III. WEAK TIME REVERSAL BREAKING 

The low energy theory relevant for the vicinity of the 
transition from an ordinary to a topological band insu- 
lator in two dimensions consists of two Dirac modes dis- 
persing around the a pair of time reversed points K and 
K' in the Brillouin zone 9 . The effective hamiltonian for 
the massive Dirac fermions near one of these points is 
given by 

H K (k) =v f k-a ± +mcr z (3) 

where <j acts on the space of doubly degenerate states at 
K, k measure the momenta relative to K and Vf is the 
Fermi velocity. The transition to the topological insula- 
tor is driven by changing the sign of the mass parameter 
m to a negative value. However we have freedom to vary 
another parameter that changes the separation between 
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FIG. 2: Numerically obtained phase diagram of the model 
in Eqs. (1,2) as a function of a staggered potential (A) and 
staggered hopping (St) in the presence of a staggered Zeeman 
field (a) or a uniform Zeeman field (b). With TRS the ordi- 
nary and topological band insulators are separated by a line 
of critical points (dashed red). When time reversal symme- 
try is broken this line can split and generically opens up into 
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a chern insulator phase of finite extent, having a xy = 
Other parts of the phase boundary can open up to become 
metallic regions of finite extent as exemplified by the gray re- 
gions in panel (b) . Notably, due to topological protection the 
gapless regions originating from the phase boundary persist 
to arbitrary large breaking of TRS. 



K and K' while keeping m = 0. In the extended param- 
eter space the transition is marked by a closed critical 
boundary surrounding the topological phase, as shown 
in Fig. (2a), dashed red line. 

It will be useful to consider the mass formally as a third 
spatial coordinate k z . Viewed in this way, Eq. (7) de- 
scribes a massless dirac cone in three dimensional space. 
The complete low energy theory including the two Dirac 
points in the extended parameter space is formally iden- 
tical to a Weyl semi- metal 10 . In contrast to the gapless 
Dirac cones in two dimensional space the Dirac points 
in the extended three dimensional space (Weyl points) 
each carry a topological index that can be formulated as 
a Chern integer 11 . This in particular protects the Dirac 
points in the extended three dimensional space from be- 
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FIG. 3: Evolution of the Weyl points in the space (k x , k y , m), 
where m is the tuning parameter, upon breaking of TRS with 
a field M. At non zero M the two Weyl points generically 
lie at different values of m. A system (II) with intermediate 
values of m is endowed with an integer Chern number from 
the two Weyl points. 



ing gapped by perturbations that break time reversal 
symmetry as long as the K and K' points are spatially 
separated. Instead the consequence of such perturbations 
is to move these points away from the plane m = (see 
Fig. 3). Moreover in absence of time reversal symmetry 
the Dirac points do not necessarily lie on the same m 
coordinate 

The physical space of a single system however corre- 
sponds to a plane of constant m. When increasing the 
value of m with time reversal symmetry broken we gener- 
ically pass two phase transitions as the plane crosses first 
one Dirac point and then the other. When the plane 
lies either above or below both Dirac points the topolog- 
ical flux emanating from them cancels and the system is 
topologically trivial. If on the other hand the physical 
plane is between the the two Dirac points the topological 
flux adds up, endowing the system with a Chern number 
of +1 or —1, depending on whether the Dirac cone with 
a positive Chern number is above or below the plane. 

The above argument implies that a small TRS break- 
ing perturbation opens the critical phase boundary into a 
Chern insulating phase. The extent of this phase initially 
grows with increasing TRS breaking field strength. At 
stronger fields the Chern insulator can eventually shrink 
and perhaps end altogether if the Dirac cones meet and 
annihilate each other. Finally we should note that the 
two Chern bands associated with the Chern insulator 
may also be separated by an indirect gap. In this case 
the Chern insulator would give way to a metal. This pos- 
sibility is realized in certain parts of the phase diagram 
plotted in Fig. (2 b). 

Before concluding this section we note that the generic 
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situation described above can be slightly modified in 
presence of various spurious symmetries. Let us give 
a few examples relevant to the model Hamiltonian (1). 
First in the case St = the operator product of time re- 
versal (T), mirror with respect to the yz— plane and s y 
is a symmetry of the Hamiltonian. As a result of this 
symmetry the two Dirac cones are constrained to remain 
on the same physical plane (m) even with broken TRS. 
For this reason the Chern insulating phases plotted in 
Fig. 2 are pinched at the points St = 0. Similarly if 
we omit the extra spin orbit term A# then the product 
of T and s z commutes with the Hamiltonian and pre- 
vents from the Chern insulator to open up at all. Fur- 
thermore, under these conditions (Xr = 0) each of the 
critical points on the line A = (at St = ±1) must have 
the the Dirac points coincide K = K' '. Then even an 
infinitesimal TRS breaking perturbation opens a gap for 
A = and therefore splits the phase boundary into two 
separated segments. 



IV. TOPOLOGICAL STABILITY OF CRITICAL 
LINES 

In the last section we discussed the perturbative sta- 
bility of the phase boundary between the topological and 
ordinary insulator to weak TRS fields. However the nu- 
merical results of section (II) show that segments of the 
critical line persist to arbitrarily strong TRS breaking. 
We shall argue that this is due to topological protection 
that is much stronger than the perturbative stability dis- 
cussed above. 

The topological invariant is rather simple to formulate 
within the framework of the non interacting band theory. 
When TRS is broken, being on a critical line in the phase 
diagram is tantamount to having a single gapless Dirac 
cone located at some point K in the Brillouin zone. If we 
go along the critical line in the two parameter phase dia- 
gram, closing a full loop parameterized by <\> G [0, 2tt], the 
Dirac point charts at the same time a closed path K(0) 
in the two dimensional Brillouin zone. We can then asso- 
ciate with the critical line, as topological invariants, the 
two winding numbers of the loop on the toroidal Bril- 
louin zone. Fig 4(a) and (c) show respectively examples 
of a topological and a non-topological loop. In the former 
case the critical line must persist in the phase diagram 
because the Brillouin-zone loop associated with it cannot 
be collapsed to a point no matter how strong the time 
reversal breaking field. 

Our goal now is to relate the winding number on the 
Brillouin zone to a more fundamental topological invari- 
ant associated with quantized pumping. This will allow 
generalization to systems with disorder and interactions 
where the single particle Bloch wave-functions are not 
eigenstates. For this purpose we consider a closed path in 
the two parameter space (A, St) that encircles the topo- 
logical insulator from the out-side, that is a path T(a) 
contained entirely in the ordinary gapped region, as de- 
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FIG. 4: The relation between windings in k— space and charge 
pumping. Panels (a) and (c) show the trajectory which the 
Weyl points make in the Brillouin zone if we move the sys- 
tem on the TI-BI phase boundary. Case (a) is an example 
with a topological winding number. Panels (b) and (d) show 
the respective trajectories of the Dirac points in the space 
(St, A, fe) for a fixed value of k\. The adiabatic path around 
the TI phase is extended to a distorted cylinder in that space. 
The charge pumped in such a cycle is given by the topologi- 
cal flux on the cylinder that originates from the Weyl points. 
This is ±2 in case (b) and zero in case (d). 

picted in Fig. (4 b). Let us now focus on the case where 
the two time reversed Dirac cones (on the phase bound- 
ary) wind in the k\ direction of the Brillouin zone. By 
fixing k\ we are dealing effectively with a one dimen- 
sional Hamiltonian, for which the pumped charge (along 
£2) is related to an integer Chern number through the 
formula 12 ' 13 

C = J2 [ ^^lm(dMd k ^ n ), (4) 

where \ip n ) = \ip n (k2,ct)) is the wave function of the ri- 
band at &2 an d and the sum is over occupied bands. 
This is an integral over the surface of the cylinder-like 
surfaces shown in Fig. 4(b,d) that measures the amount 
of topological flux emanating from the Dirac points en- 
closed by it. These Dirac points can only reside on the 
TI-BI phase-boundary. But for a given value of k\ we are 
not in general guaranteed to have any. For example, if 
the Dirac cones at K and K' chart paths in /c-space that 
do not wind on the Brillouin zone (see e.g. Fig 4(c)), 
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then we can fix k\ to a value such that the effective one 
dimensional Hamiltonian along k^ is fully gapped, even 
on the phase boundary. On the other hand, if the Dirac 
cones at K and K' each winds around the Brillouin zone 
in the direction of k\ then for any fixed value of k\ we 
must have one Dirac point on the phase boundary origi- 
nating from the trajectory of K and the other from K' . 

The Dirac points originating from K and from K' 
contribute equal rather than opposite flux to the sur- 
face integral, contrary to what might have been naively 
guessed. To see this first recall that each Dirac point 
that is enclosed by our surface can be viewed as a hedge- 
hog singularity in the relevant three dimensional space 
k = (&2,££, A), with low energy hamiltonian near the 
Dirac point i^K,fci(^) = —Sk-a. The contribution of the 
Dirac point to the flux through the surface is the integer 




where k = Sk/\Sk\. This is an integer that must re- 
main fixed as a function of the continuous parameter k±. 
Let us then compare the flux Ck,^ emanating from the 
Dirac point at K to that from its time reversed partner 
Ck' -fc r Since in the Hamiltonian —5k • a both k\ and 
cr y , and only those, change sign under time reversal sym- 
metry the topological flux emanating from the two time 
reversed points is equal. Studying the winding number 
of the phase boundary between the TI and BI, we found 
it to be ±1 along the direction of the staggered hopping. 
Therefore adiabatically changing the parameters of the 
model (1) to encircle the topological insulating phase, 
while staying in the ordinary insulator, entails pumping 
of exactly two electrons per row of the lattice. Such quan- 
tized pumping has been shown to be robust to disorder 
and interactions that invalidate the simple band picture 
13 . We therefore expect the gapless phase boundaries to 
persist even in presence of interactions . 

V. SUMMARY AND DISCUSSION 

We showed that the critical phase boundaries separat- 
ing a two dimensional topological insulator from an or- 
dinary insulator are more robust to TRS breaking than 
the bulk phases themselves. First we established the per- 
turbative stability of the critical lines in the generic case 
through a mapping to a Weyl metal in one higher dimen- 
sion. A corollary of this analysis is that the critical phase 
boundary opens up into a Chern insulating phase or into 
a metallic phase with finite extent in the phase diagram. 
Next we formulated a criterion for yet stronger topolog- 
ical protection of the gapless phase boundaries that is 
insensitive to interactions and disorder. When the topo- 
logical criterion holds the phase boundaries can split into 
segments, but the segments must persist to arbitrarily 
strong TRS breaking fields. 

It is interesting to point out a surprising connection 
between the observations made here and previous work 



on the ionic Hubbard model (i.e. Hubbard model sup- 
plemented by a staggered potential). The latter is used 
as a model to describe transitions between a band in- 
sulator and a correlated Mott insulator with increasing 
interaction strength. The onset of the Mott phase is typ- 
ically defined through the emergence of magnetic order, 
both phases are however gapped to single particle exci- 
tations. A surprising and yet unexplained observation 
made in a number of numerical studies is that the mag- 
netically ordered state itself is split into two regions sep- 
arated by a metallic region with vanishing single particle 
gap 14 16 . Using a simple Hartree Fock description, the 
metallic phase can be seen to be smoothly connected to 
a phase boundary between a TI and ordinary band insu- 
lator at zero magnetization and non-vanishing spin orbit 
coupling (see App. A). The persistence of the metal- 
lic phase deep into the magnetically ordered state, even 
without spin orbit, is due to the topological robustness 
of the phase boundary. 

Finally, as an extension to this work it would be in- 
teresting to classify the robustness of phase boundaries 
between topological and ordinary insulators in higher 
dimensions (as well as higher dimensional parameter 
spaces). This can be done using the general classifica- 
tion of defects and pumps in such phases 17 ' 18 . 

Appendix A: Protected gapless regions in the ionic 
Hubbard model 

In this appendix we relate the metalllic regions that ap- 
pear inside the magnetically ordered phase of the ionic 
Hubbard model 14 16 > 19 with the persistent phase bound- 
aries discussed in the main text of this paper. The ionic 
Hubbard model is defined by the following Hamiltonian 
on a square lattice 

H e l = tC L C j,s + A & n M ( A1 ) 

(ij),s i,s 

+ ^C/(n a -l/2)(n a -l/2), 

i 

where Ci jS annihilates an electron with spin s on the i'th 
site, & = ±1 denotes an alternating sign on the two 
sublattices and n^ s is the occupancy. 

The qualitative features of this model can be seen from 
a Hartree Fock (HF) analysis. The interaction term is 
decoupled assuming a staggered magnetic moment 

Y^U(n itt - 1/2) (n u - 1/2) « M^&ctX,,'*,*', 

i i 

(A2) 

and the resulting quadratic Hamiltonian is solved to- 
gether with the self consistency condition 20 

Mi = ^(4,^,s /C j>')m- (A3) 
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FIG. 5: A square lattice can be viewed as a deformed hexag- 
onal lattice with one extra bond per side (dashed line). 

The mean field theory shows a second order transition 
into the (AF) magnetically ordered state at a critical 
interaction strength U c oc t/ln(t/A). The single parti- 
cle excitation spectrum of the mean field Hamiltonian is 
given by 

E Ks = ± ] /e 2 k + (sM stg + AY, (A4) 
€k = 2t(cos(k x a) + cos{k y a)), 

where s denotes the spin along the x— direction. Notably 
as the magnetization increases it eventually reaches the 
value \M stg \ = |A|. At this point the excitation gap for 
spins with s = —sign(A/M stg ) closes and a gapless point 
appears. 

We argue that the gapless point inside the magnetically 
ordered phase can be related to the critical point between 
a TI and BI in presence of TRS. To this end we add the 
following perturbation which gradually deforms the free 
part of ionic Hubbard model on the square lattice into 
the Kane and Mele model defined on a honeycomb lattice 

21 

V = i\ V i3 S As z C 3s z , (A5) 

- x c i z c js z . 

(ij)v,s z 

Here ((ij)) denotes pairs of sites connected by two near- 
est neighbor bonds which do not include the dashed-bond 
in Fig. (5) and (ij) v denotes sites sharing a dashed- 
bond. The spectrum of the mean-field Hamiltonian to- 
gether with the addition V is given by 

E kt8 = ±Vl^| 2 + M + A)2, (A6) 
e k = 2te ik y a (cos(k x a) + cos(k y a)) - Ae 2i ^ a , 

A = x U\ 2 [sin(2k x a) - 2cos(k y a)sin(k x a)] 2 + M 2 tg , 



where s = ±1 now denotes the spin projection along the 
direction (M stg , 0, 2X[sin(2k x a) — 2cos(k y a)sin(k x a)]). 



035 



A 




0.5-1 1.5 

M x 



FIG. 6: Regions with gapless charge excitations in the ionic 
Hubbard model at t = A = 1 as a function of the sponta- 
neous magnetization M s t g and the spin-orbit perturbation A 
(see Eqs. (A1,A2,A5)). The A = axis corresponds to the 
usual ionic Hubbard model where a gapless point appears for 
M stg = A. Increasing A at M stg = another gapless point 
appears marking the transition between a BI phase and a TI 
phase. Interestingly, these two points are smoothly connected 
implying that the former one could be seen as part of the pro- 
tected gapless surface discussed in the main text. 



Using Eq. (A6) we obtain the equations for the criti- 
cal points where the gap closes: 



4X 2 [sin(2k x a) - 2sin(k x a)} 2 = A 2 - M 2 ^, (A7) 
cos(k x a) = (^- - 1). 



These equations give a critical line in the space A versus 
M stg which is shown in Fig. 6. The critical line extends 
all the way to the axes at A = corresponding to the 
ionic Hubbard model on the square lattice. Thus the 
gapless point in the magnetically ordered phase of the 
ionic Hubbard model can be understood as persistence 
of the critical surface separating the TI and BI phases 
due to the topological protection discussed in the main 
text. 
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